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Abstract. Solution of the Navier-Stokes equations with initial conditions (Cauchy problem) for 2D and 3D cases was 
obtained in the convergence series form by iterative method using Fourier and Laplace transforms in paper T. For several 
combinations of problem parameters numerical results were obtained and presented as graphs. 

Estimated formula for the border of the parameter area of convergence of the iterative method was obtained in paper [2] . 

This paper describes numerical proof of convergence of the iterative method for solution of the Cauchy problem for the 
Navier - Stokes equations. Usage of estimated formula for the border of the parameter area of convergence of the iterative 
method is shown for wide ranges of the problem's parameters. 

1. Introduction. 

Example of the solution of the Navier-Stokes equations with initial conditions (Cauchy problem) with a 
particular radial applied force for 2D case was obtained in paper [T]. 
We have initial conditions 

(1.1) u{x,0) = u^{x) = (xeR^) 

and a radial applied force 



(1.2) /i,(f,t,T,^) = F„r»+ie 



(1.3) /i^(f,(^,r) = 



Here {f,ip) are the polar coordinates; n is a circumferential mode, n=l,2,3,...; Fm^n - constants; v is & 
positive coefficient of the viscosity; t, t - time. 

A proof of convergence of the iterative method for solution of the Cauchy problem for the 3D Navier - 
Stokes equations was described in paper [2j . As one of the results of this proof an estimated formula for the 
border of the parameter area of convergence of the iterative method was also obtained: 



It is easy to show for an estimate solution of the Cauchy problem for the 2D Navier - Stokes equations 
that estimated formulas for the border of the parameter area of convergence of the iterative method for 
different circumferential modes have the same form: 



(1.5) 4^<1 n = 1,2,3, 
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2. Numerical analysis based on estimated formula. 

Numerical results for circumferential modes n = 1, 2, 3, 4, 5 with a radial applied force (|1.2p are shown 
on FIG. 1 - 5. 

We plot curves of the form: 



in coordinates (horizontal, logarithmic scale) and /i (vertical) for various v. 

Plots for V = 0.01, 0.1, 0.3, 0.75, 1.0, 1.5 are displayed on FIG. l.a, 2.a, 3.a, 4.a, 5.a. 

Areas of convergence of iterative method are above these curves. 

For each circumferential mode n results are calculated for amplitudes of Fn taken by formula 



(2.2) Fn = 10''/n A: = 0,1, 2, 3. 

Amplitudes of F„ are selected from range 0.2 < F„ < 1000. 

Since areas of convergence of the iterative method are above displayed curves, we have selected the 
highest curve from our set of tests for further analysis. For all circumferential modes n and all chosen values 
of Fn the corresponding values of fj, are appeared to be in a range 2.0 < /i < 18.0. All of them are shown as 
black dots on the curve for ly — 0.01. 

For selected values of F„ , fj, and comparison of calculated Maximums of Velocity on the first step | ui \ 
(red plots) and the second step | | (blue plots) of iterative process is displayed on FIG l.b, 2.b, 3.b, 4.b, 
5.b. The upper graph is shown in coordinates Fn (horizontal, logarithmic scale) and Amplitude of Velocity 
(vertical, logarithmic scale) for ly — 0.01. The bottom four graphs are plotted in coordinates v (horizontal) 
and Amplitude of Velocity (vertical, logarithmic scale). For all graphs we see that | | << \ ui |. On 
FIG l.c, 2.C, 3.C, 4.C, 5.c the comparison of first step amplitudes and second step amplitudes is displayed 
in plane (p — [0, tt] for various radii r in a range < r < 3. Also from these graphs we see that for all 
circumferential modes amplitudes of velocity on the second step of iterative process are smaller than corre- 
sponding amplitudes on the first step. Left graphs have full amplitude scale and right graphs are zoomed 
into low amplitude values for better display of second step amplitudes. More detailed results of amplitude 
comparison are shown in [l] . 

From all the graphs below we may conclude that the iterative process described in [T] and [5] has very 
good convergence. For higher values of fi with the same values of Fn the convergence of the method will be 
even better. 
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FIG.4.a. 
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